Gluon fusion production of the Higgs boson 
in a calculable model with one extra dimension 



Giacomo Cacciapaglia a '[], Marco Cirelli a, 0, Giampaolo Cristadoro 6 '^ 

a Scuola Normale Superiore and INFN, Piazza dei Cavalieri 7, 1-56126 Pisa, Italy 
b 'Enrico Fermi' Department of Physics, University of Pisa, 
via Buonarroti 2, 1-56126 Pisa, Italy 



Abstract 

In an extension of the Standard Model with one extra dimension and N=l supcrsymmetry com- 
pactified on R 1 /Z2 x U 2 , we compute the Higgs boson decay width into two gluons, relevant to 
Higgs production in hadronic collisions. At one loop, the decay width is significantly suppressed 
with respect to the SM. For a compactification radius R = (370 ± 70 GeV) -1 and a Higgs mass 
mjj = 127 ± 8GeV, as expected in the case of a radiatively generated Fayet-Iliopoulos term, we find 
it to be less than 15% of the SM result. 
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1 Introduction 



Recently, a large interest has arisen about the idea that compact extra dimensions can be not so small 
as previously thought Q. The main drawback, from a phenomenological point of view, is that most of 
the models based on this idea lack quantitative aspects, so that it is often hard to find a connection with 
phenomena that take place at a known physical scale. The model we consider [0], instead, leads to the 
calculability of several observables, due to the underlying supersymmetric structure of the theory. In 
this framework, the Higgs mass fU, the branching ratio B — > X s ^ || and the muon anomalous magnetic 
moment jjj have already been computed and found insensitive to the UV physics. 

Following this line, in this paper we compute the Higgs decay width into two gluons, with the aim of 
presenting a reliable and cut-off independent calculation in extra dimensions. The main interest about 
this process is that it can be related to the cross section for Higgs production |j] through gluon fusion. 
At a large hadron collider (LHC), this is expected to be the most important channel for light SM Higgs 
[||. Indeed, in the model we consider, a light Higgs is predicted in the range ran ~ 120 -j- 160 GeV . 
Moreover, we focus only on zero mode gluon fusion since we expect the parton distribution of KK gluons 
in the proton to be negligibly small. 



2 The model 

The framework of our calculation is based on the model of ref . Q . This has the same matter content 
and gauge structure of the Standard Model, but is embedded in a 5-dimensional space time with 
N=l supersymmetry. The orbifold compactification of the extra dimension on R 1 /Z2 x Z 2 provides 
supersymmetry breaking a la Scherk-Schwarz |7]]. For our purposes the compactification radius R is a 
parameter and we neglect any dynamical mechanism that stabilizes it. The presence of only one Higgs 
hypermultiplet constrains the form of the scalar potential. As pointed out in ref. a Fayet-Iliopoulos 
(FI) term £D can be located on the two branes, thus introducing a free parameter in the potential. 
Nevertheless, after ElectroWeak Symmetry Breaking driven by radiative corrections, the Higgs mass is 
predicted to lie in the range 120-^160Gey, while the value of 1/R, which is correlated to the Higgs mass, 
can increase up to about a TeV ||. In the particular case that the FI term is radiatively generated, 
with a cut-off at a scale ~ 5/R, its effect is small enough and the predictions of ref. [|| still hold: 



m H = 127 ±8 GeV 
l/R = 370 ± 70 GeV. 



In [10] a hypercharge current anomaly is found for the model under consideration. It is localized on the 
two branes and has the property that the integrated anomaly vanishes. Gauge invariance can however 
be recovered if the theory is modified in a suitable way. The overall consistency of the model with the 
addition of a Chern-Simons term, for instance, is currently under examination. Anyhow, we believe 
that our calculation should not be affected by the fixing mechanism. At one loop, indeed, the Higgs 
decay is only determined by the brane interactions and does not depend on the (bulk) gauge electroweak 
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structure. For this reason, also, our results can be extended to a broader class of models: for instance, 
in a model with two Higgs hypermultiplets [11] the same calculation holds after a proper identification 
of the Higgs field involved in the vertex. 



The relevant lagrangian for our calculation comes from the quark Yukawa interactions, located on 
the branes. We remind that the vertices of gluon emission are the same for all the KK states, with 4- 
dimensional coupling g s = / ^/2ttR. In terms of mass eigenstates, for every flavour q the Lagrangian 
is: 

oo 

c int = - E { « OM 2 H^n + $ n j± + Kom r^t + h-c) } 

n=l 

- (mo(fc) V^o + h-c) (1) 

where mo, rrv^ and m^ are the field dependent masses for quark and squark towers. They have the 
same form for every flavor q, even if the mass matrices arising from Yukawa couplings, y q , located on 
different branes are not the same (see App. [A]): 



in 



± 2n-l 



n 

± _ ™ _i_ 2n 

n 



m q (<p H ) ± f (2) 
m = m q ((j) H ) = ^arctan (f y q R<f} H ) 

Substituting for the complex field 4>H = v + H/\^2, this provides us with the full expression of the 
couplings of the Higgs field H to the mass eigenstates. Taking the linear terms in H and keeping all 
orders in m q R = e q , we find: 



oo oo 

Cint = - A, %r H -\ q J2 ^t c H + h.c. - \ q (±)2m» ($±t$± + $f*$^)H 

n=l n=l 



(3) 



where 



J_ 1 [Gf 2V2 tan(f eg) 

9 ^l+tan 2 ^) V ^ l + tan 2 (fe g ) U 

Note that the Yukawa coupling is related to the Fermi constant -7= = ^ in a non trivial way, due to 
the mixing with the massive KK states. The standard expression, X q = ,J^=y/2m q , is recovered for 
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3 The calculation 

The graphs that contribute at one loop are shown in figure |]. For a color triplet scalar </> and a color 
triplet fermion ^ with mass m s and mj respectively, coupled to the Higgs field via C = Xf^f^fH + 
\ s (j)'(j)H, the amplitude can be written as: 
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Figure 1: Feynman diagrams contributing at 1-loop to H — > gg. Dashed lines stand for scalars, 
full lines for fermions. 



A ab = -^ T [(Pi • ^(Pa • ei) - (pi • P2)(ei • e 2 )] ^ 
where £j and j>j are the gluon polarizations and momenta, and 

-1 rl-x 



(5) 
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Let us first review the Standard Model result [JT^]. In the present context it comes from the zero 
mode fermionic graph for Xf = y q , and it amounts to: 
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6 v ^-^ 



(7) 



where £, 



2 2 

p m H 

q o o": 



and the function g is 



A 5A /(0 = ^{2£+(£-4) 



9(6 



2arcsin 2 ( 



log 



n-yi-^- 1 



i7T 



£ < 4 

e > 4 



The decay width is related to the amplitude, after summing over polarizations and color indices of the 
final state gluons. In the SM it is: 
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Fsm(H -> gg) 



Gf oi s m 6 H 
^2 36vr 2 



(8) 



The main contribution comes from top loops, while the other quarks contribute to a few percent. 
The dependence on the Higgs mass, apart from the m? H factor in To, is encoded in £, in such a way that 
the functions Asm run from « 0.9, for mn = 120 GeV , up to ps 1 for m/r = 160 GeF. 



Let us now move to our model. The couplings are given in eq. (0), and we consider contributions 



from top and stop KK towers. Defining e = m t R and £ = — ^, the result can be written as: 



HIT 



T f = 8\ t R J dxdy (1 - Axy) { ^J^ + V, ^ e ) } 
JF S = 4Afi? y dxdy 4xy s(x, y, £, e) x 2 



(9) 
(10) 



where the factor two in ,F S takes in account the two KK towers of scalar partners and the functions / 
and s contain the contributions of the KK modes: 



/ 



e + 2n 



2n 



^ {(2n + e) 2 - £e 2 xy + (2n - e) 2 - £e 2 xy 
^ / e + 2n - 1 e - 2n + 1 



■J \ (2n — 1 + e) 2 — £e 2 xy (2n — 1 — e) 2 — ^e 2 xy 



(11) 



Let us remind that the Higgs mass depends on R, so that in our model the two parameters e and £ are 
related to each other. 



3.1 £ -> limit 

To compute the KK state contribution exactly is not a trivial task. As a first step we note that the 
sums in / and s simplify in the limit £ — > 0. In this case we are neglecting contributions of relative 
order 0{mnR) 2 ■ These contributions may be important expecially for large R, where ran ~ 120 GeV 
and then itlhR > 0.4. However, we use this limit just to have a reliable estimate for small radius 
(1/R ~ 1 TeV), where ran ~ 160 GeV and then mnR ~ 0.15, and to have an analytic control on the 
calculation. 

In this case, the result is: 



T(H - gg) = T 



tan 



vre 1 + tan 2 (fe) 



A t - 1 + ^7recot 



1 /7T 

-7retan — e 
4 V2 



+ A 



(12) 
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where A t>b = A S M{Ct,b)- 

As we can see from figure the dependence on the Higgs mass in eq. (|l2|) is smooth in the range 
120 -7- 160 GeV and then the correlation between ran and R can be consistently neglected. The decay 
width is significantly different from the SM one and, in particular, it vanishes for 1/R ~ 300 GeV (that 
corresponds to ran ~ 120 GeV ||). This means that the contributions of the top and stop KK modes 
cancel the SM result. It is relevant that this happens for every value of A b ^\ 

Also notice that for the unphysical value 1/R = nit the decay width diverges. This is a remnant of 
the approximation we are working in, as corrections C(m^i?) 2 are relevant in this region. 

3.2 Full calculation 

In the previous paragraph, the dependence on ran has been kept only in the SM result. This is, as we 
saw, an approximation less appropriate for small values of 1/R. In order to evaluate how the result 
fully depends on ran, one could keep the £ dependence in eq. (|ll]) and then expand in the parameter e. 
But, two issues arise. The first one is that if we expand eq. ([l2|), the series is oscillatory, and stabilizes 
only at very high order in e. The second point is that for small 1/R the parameter e is too large and 
the expansion is not meaningfull. 

What one can rather do is to perform the sums in / and s and then evaluate numerically the residual 
integrals on Feynman parameters x and y. The functions / and s can be simplified to: 

1 We have checked that working out a formula similar to eq. ( |l2| ) including also bottom KK towers, there is no significant 
change since the parameter th is quite small. 
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where tp(z) 



2+ £ (v^-l 



2+e(V^+l) 



-2eY I = 

n=l v v ' 

oo 1 

-2eV L_ = 

^ (2n - 1 + eV^) 2 - e 2 

r'(2;)/r(z) is the polygamma function of order zero. Now, the net result is 



* ( 1+ <^- X ) 



^ ( 1 + <V^+1) 



(13) 



r r 



tan(fe) 



vre 1 +tan 2 (fe) 



A t + 3eE(£,e) 



+ A, 



where 



1 rl—x 

dx / [(1 — 4xy)f + 4xys] . 

./o 



(14) 



Eq. ( |l~4|) contains the exact dependence on the related parameters e and £. In order to estimate 
how the approximation in the previous paragraph works, we fix m# = 127 GeV and in fig. || we plot 
the result as a function of R compared to eq. ( |l2| ) in the same limit. The two curves differ only slightly 



in the region 1/R > 300 GeV. Moreover, for small values of 1/R the apparent divergence of eq. (^L2|) is 
not present any more. 

In figure |3], we can also read the value of the decay width for 1/R = 370 ± 70 GeV. To do this, 
we first have to take into account the effect of the brane z factors discussed in ||, that are the main 
source of error on 1/R. Eqs. (p^) and (|l4l) are actually functions of R = R/(l — z) instead of R, since 
the dependence on the radius only comes from the top and stop mass eigenvalues. It is R which is 
determined by the minimization of the Higgs potential, so that its uncertainty does not come from the z 
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factors but from other sources. It is therefore estimated to be smaller, at 10% level ||. Thus, it follows 
that the decay width, in the case of a radiatively generated FI term, is at 5% -j- 15% level of the SM 
one and is only affected at quadratic order in z. 



4 Conclusions 

In this paper we have calculated the full one loop amplitude for the decay width of the Higgs boson 
into two gluons for any value of the compactification radius. We find a result significantly suppressed 
with respect to the SM value. This is due to a partial cancellation between the SM and the top and 
stop KK tower contributions. In particular, this cancellation is complete and the cross section goes to 
zero for a value of 1/R « 300 GeV. 

For the value of the compactification radius of 1/R ~ 370 ± 70 GeV, as expected in the case of 
a radiatively generated FI term, the one loop cross section is reduced below the 15% level of the SM 
prediction. Notice, however, that the two loop amplitude, of order a 2 , is expected to be relevant, since 
in the SM [jl3| it is as large as 40% of the 1-loop one. 



A Mass matrices 

In ref . |J , the diagonalization of the mass matrices arising from Yukawa couplings located on the brane 
y = is presented. Dealing with couplings located on the brane y = some subtleties appear. First 
of all, the matrices are, for scalar s and fermions: 



M 2 + 4m 2 r]E 
-2mM ■ E 



{^T,0,^T,l,1pu,l) 



( 



m 







—2mE ■ M \ 




M 2 ) 


\ &T,k 


V2mJ T \ 




2mJJ T M 


lpU,k 


M j 




Eik = {-l) l+k and m = 



(15) 
(16) 

where r, = Y, k 1, M lk = |<5 Zfc , M lk = '^5 lk , J k = (-l) fe , E lk = (-l) l+k and m = y q v. The difference 
with ref. || is the presence of some —1 signs that come from the values of the wave functions of the 
fields, whose dependence on y in given by their parity. But, carefully performing the calculation, it is 
clear that the eigenvalues coincide with those produced by the matrices on the brane y = 0, eq. 
Note also that in this case <f> c is the field in the N=l supermultiplet of ip. The effect of these signs 
is physically relevant only when the two branes interfere. This happens when the observable receives 
contributions from both the branes (see for example 0). 

The calculation of these matrices deserves another comment. In fact, if we write the bulk lagrangian 
as in W% and add the two brane interactions as in H, it is not trivial how to get rid of the auxiliary 
fields F and F c . For F it is straightforward since both the bulk and the y = brane lagrangians are 



written in terms of the same multiplets, while for F c we refer the reader to [15]. The result is that the 
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term 5(y — TTR/2)d^4>-^W' ', where W is the superpotential localized on the y = irR/2 brane, has a 
minus sign with respect to the analogous term generated by F. 
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